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Electronic and magnetic properties of Gai-a^Mn^^As, obtained from first-principles calculations 
employing the hybrid HSE06 functional, are presented for x = 6.25% and 12.5% under pressures 
ranging from 0 to 15 GPa. In agreement with photoemission experiments at ambient pressure, we 
find for x = 6.25% that non-hybridized Mn-3d levels and Mn-induced states reside about 5 and 
0.4 eV below the Fermi energy, respectively. For elevated pressures, the Mn-3d levels, Mn-induced 
states, and the Fermi level shift towards higher energies, however, the position of the Mn-induced 
states relative to the Fermi energy remains constant due to hybridization of the Mn-3d levels with the 
valence As-4p orbitals. We also evaluate, employing Monte Carlo simulations, the Curie temperature 
(To). At zero pressure, we obtain Tc = 181 K, whereas the pressure-induced changes in Tc are 
dTc/dp = 4-4.3 K/GPa for x = 12.5% and an estimated value of dTc/dp « 4-2.2 K/GPa for 
X = 6.25% under pressures up to 6 GPa. The determined values of dTc/dp compare favorably with 
dTc/dp = +(2-3) K/GPa at p < 1.2 GPa found experimentally and estimated within the p-d Zener 
model for Gao.gsMno.oTAs in the regime where hole localization effects are of minor importance [M. 
Gryglas-Borysiewicz et al., Phys. Rev. B 82, 153204 (2010)]. 


I. INTRODUCTION 

The application of hydrostatic pressure serves as a 
unique tool for tuning the various essential parameters 
in semiconductors [l|. High-pressure experiments and 
theoretical studies were important, for instance, in the 
delineation of the band structure of kw semiconductors 
like Si, Ge, and the III-V compounds [3. This approach 
was also broadly used (both in theory and experiment) in 
the studies of bistability of donors {DX centers) in GaAs 
and GaN [|. 

For spintronic semiconductors, high-pressure measure¬ 
ments can play an important role in understanding the 
mechanisms that are responsible for the coupling be¬ 
tween magnetic ions in dilute magnetic semiconductors 
(DMSs): the changes in volume of the solid influence 
the local exchange interaction between magnetic dopants 
providing information about the mechanism responsible 
for the observed macroscopic magnetization Q. Pres¬ 
sure can also help to elucidate the relationship between 
the coupling mechanism between magnetic dopants and 
the underlying band structure of the DMS. The states in¬ 
troduced by the dopants may have characteristics of the 
band edge states and follow the band edge directly when 
the pressure is applied or may tend to be highly localized, 
therefore, influenced by the entire Brillouin zone (the 
pressure dependence of deep levels does not follow any 
particular band edge) Bearing in mind the sensitiv¬ 
ity of the exchange interaction between magnetic dopants 
to interatomic distances in DMSs, their experimental and 
ab initio theoretical studies under applied pressures are 
the most direct way to test different theoretical models 
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that have been proposed to emlain, for instance, the fer- 
romagnetism in (Ga,Mn)As Q. 

There exists a number of experimental studies on the 
variation of Tq under pressure for DMSs: (Pb,Sn,Mn)Te 
0, (In,Mn)Sb M, (Sb,V)2Te3 0, (Ga,Mn)As 
and (Ge,Mn)Te Ij. The results of those studies reveal 
that for samples with high concentrations of band holes, 
Tc increases with an increasing hydrostatic pressure p 
according to the expectations of the p-d Zener model 
[1^ . However, there are two worthwhile exceptions, i.e., 
the systems showing dTc/dp < 0. The first is a narrow- 
gap topological insulator (Sb,V)2Te3 0, in which rather 
than intraband excitations (i.e., the standard Ruderman- 
Kittel-Kasuya-Yosida mechanism), inter band excitations 
(i.e., the Bloembergen-Rowland interactions) give a dom¬ 
inant contribution to the comling between localized spins 
of transition metals (TMs) 0. The second is (Ga,Mn)As 
in the vicinity of the metal-insulator transition, where the 
increase in the p-d hybridization enhances hole localiza¬ 
tion and, hence, diminishes Tq 0- 

Not much has been accomplished from the ab ini¬ 
tio side so far. The pressure dependence of Tc for 
Gao.gsMno.osAs was studied by Bergqvist et al. [l^ 
within the local-density approximation (LDA) and the 
LDA+ U approach for pressures slightly crossing 7 GPa 
for which the direct-to-indirect band gap transition oc¬ 
curs for GaAs within the LDA 0- It was established 
(employing the LDA+R) that Tq increases for the de¬ 
creasing lattice constant of Gag .gsMno, 05AS. This behav¬ 
ior was more pronounced within the mean-field approx¬ 
imation (MFA) than for the Monte Garlo (MG) simula¬ 
tions. 

Here, we present first principles calculations for 
the structural, electronic, and magnetic properties of 
Mna;Gai_a;As with x = 6.25% and 12.5%, exposed to 
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pressures from 0 to 15 GPa. Therefore, we cover the 
whole range of pressures for which the zinc-blende struc¬ 
ture of GaAs is stable M- Although the focus of this 
paper is to study the pressure dependence of Tq, some 
other issues are also discussed, such as the influence of 
pressure onto the electronic and magnetic properties of 
Gao.9375Mno.o625As, the energetics of the formation of 
Mn pairs in GaAs, and finally, the pressure dependence 
of the exchange coupling between Mn ions both at purely 
substitutional and mixed substitutional-interstitial sites. 

The paper is organized as follows. In Sec. II the com¬ 
putational approach is discussed. In Sec. Ill the results 
are presented, beginning with a discussion in Sec. IIIA 
of the structural, magnetic, and electronic properties of 
the studied DMS under pressure. The results concerning 
the cluster formation and the pressure dependence of Tq 
are presented in Secs. IIIB and IIIG, respectively. We 
end with some concluding remarks in Sec. IV. 


II. COMPUTATIONAL DETAILS 

The first-principles spin-polarized calculations are per¬ 
formed using a plane-wave basis and Troullier-Martins 
norm-conserving pseudo-potentials (l^ . [T^ as imple¬ 
mented in the Quantum-ESPRESSO package [l^. We 
employ the hybrid HSE06 exchange-correlation func¬ 
tional [l^ [23 that appears particularly suitable for sys¬ 
tems containing highly correlated localized electrons on 
TM d shells and itinerant band carriers. The plane-wave 
cutoff is set to be 40 Ry and a 8 x 8 x 8 fc-point mesh 
for the Brillouin zone sampling is used. The calcula¬ 
tions are done using a body centered cubic lattice of su¬ 
percells containing 32 atoms. Its primitive vectors are 
ai = a (—1,1,1), 3,2 = a (1, —1,1), and = a (1,1, —I), 
where a is the edge of the conventional unit cell. By re¬ 
placing one Ga ion with a Mn ion, we obtain a ferromag¬ 
netic (EM) DMS with the Mn concentration of 6.25%. 
We also study a system with two Mn ions in the Ga sub¬ 
lattice (it corresponds to Mn concentration of 12.5%). 
Finally, for completeness, we study a system with one 
substitutional and one interstitial Mn ion in a chosen su¬ 
percell. The calculated equilibrium values of lattice con¬ 
stant, bulk modulus, and pressure derivative of the bulk 
modulus for GaAs are a = 5.70 A, R = 69.5 GPa, and 
B' = 4, respectively, and agrees well with earlier theo¬ 
retical and experimental reports [2l|, [2^ . Moreover, the 
computed GaAs lattice constant is close to the experi¬ 
mental values, 5.68 A (5.70 A) for Gai-^Mn^jAs with 
X ss 6% {x ~ 12%) [23. The Murnaghan equation of 
state is used to estimate the theoretical pressures. Since 
HSE06 calculations are computationally very demanding 
in comparison to PBE ones, the internal coordinates of 
atoms of GaAs doped with Mn were optimized at the 
PBE level of theory. 

Structural optimization is a prerequisite to a successful 
computation of the electronic structure. This is because 
relaxation around the Mn impurity could, in principle, 
influence the electronic structure substantially. We have 


explicitly allowed for structural relaxation by displacing 
each atom from its equilibrium position for unit cells at 
zero and higher pressures. We find, for instance, that at 
zero pressure, the introduction of one Mn ion has a large 
effect on the position of the first shell of As neighbors, 
where the Mn-As bond length is elongated by 3.5% with 
respect to the Ga-As bond length of 2.469 A in GaAs. 
The tetrahedral symmetry is preserved at all pressures. 

The MG simulations are performed using the VAM¬ 
PIRE software package (IJ. The values of Tq are cal¬ 
culated from the thermal average of magnetization vs. 
temperature curve. For these calculations, the cell size 
is chosen to be 22 x 22 x 22 and for each temperature, 
we perform 240000 equilibration steps and the averag¬ 
ing loop is set to 240000 steps, which gives reasonable 
averaging. 


III. RESULTS AND DISCUSSION 
A. Structural, magnetic, and electronic properties 

For Gai_a;Mna;As with x = 6.25% and zero external 
pressure, we find a total magnetic moment of 4.27 Bohr- 
magnetons (/is) per super-cell. This value is reduced 
comparing to the magnetic moment of 5/rB for the free 
Mn atom by spins of band holes which are aligned accord¬ 
ing to the antiferromagnetic p-d exchange interaction. 
We note that the orbital contribution to the magnetic 
moment, brought about by the spin-orbit interaction, is 
below 0.1/iB per Mn ion in (Ga,Mn)As For larger 

pressures, the total magnetic moment does not change 
significantly, for instance, for 7 and 15 GPa, it assumes 
the values of 4.23 and 4.19/rB, respectively. 

Galculations of the density of states are presented for 
hydrostatic pressures of 0, 7, and 15 GPa in Fig. [TJ As 
seen, the majority and minority spin components show 
a band-gap, indicating that the introduction of the Mn 
ions does not destroy the semiconducting nature of the 
material. This is further shown in Fig. [5] where we plot 
the pressure dependence of the majority spin band gap 
compared with that of pure GaAs. In the isolated im¬ 
purity limit, Mn in GaAs is a shallow acceptor, situated 
~0.1 eV above the valence band maximum [^. It hy¬ 
bridizes primarily with the valence band, and completely 
merges with it for a few percent of Mn (no impurity 
band). From Fig. (TJ it is apparent that the introduction 
of Mn to GaAs results in the Mn-induced spin-splitting 
of the valence band, where the valence band maximum of 
the majority spin is ^0.3 eV above that of the minority 
spin. This prediction is in qualitative but not quanti¬ 
tative agreement with experiment [2^, since experimen¬ 
tally a smaller splitting of ~0.2 eV is expected for x ss 6% 
0. It is, however, documented in the literature that the 
existing implementations of the density functional theory 
|I^T) tend to overestimate the p-d exchange splitting 

From Fig.lll we can also see that the Fermi level is not 
located entirely within the majority-band. This means 
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FIG. 1. (color online). Spin-resolved DOS in logarithmic 
scale of Mn 0 . 0625 Ga 0 . 9375 As (one Mn atom per supercell) at 
hydrostatic pressures of 0, 7, and 15 GPa. Light gray color 
denotes the total DOS whereas the Mn-3d projected DOS is 
shown in black color. All the energy values are referred to the 
position of the main peak of the Mn-3d states. The dashed 
vertical line indicates the Fermi level. 


that the system is not completely half metallic, and as a 
consequence, the total magnetic moment per Mn ion is 
larger than 4/iB. However, since our approach overesti¬ 
mates the exchange splitting, the real magnetic moment 
is closer to S/tb than our theory predicts. The Mn-3d 
states in Fig. [1] are for zero pressure located 5.1 eV be¬ 
low the Fermi level. This value is slightly overestimated 
with respect to recent experimental results in which the 
main peak of the Mii-3d states was observed at 4.5 eV 
below the Fermi level [13 ■ However, it should be noted 
that the theoretical results for the location of the Mn- 
3 d states strongly depend on the employed theoretical 
approach. The present results significantly improve pre¬ 
vious DFT calculations that used the standard LDA or 
general gradient approximation (GGA) [llj. In addition 
to the Mn-3d states, we also find, in agreement with ex¬ 
periment [13, Mn-induced states centered about 0.4 eV 
below the Fermi energy. For elevated pressures (7 and 
15 GPa in Fig. [1]), we obtain a pressure-induced shift 
of the Mn-3d levels, Mn-induced states, and the Fermi 
energy towards higher energies, whereas the position of 
the Mn-induced states relative to the Fermi energy re¬ 
mains constant. This is an additional indication for hy¬ 
bridization of the Mn-3d levels with the valence As-4p 
orbitals. It is also clear from our computations that the 
Mn-derived spin-polarized feature in the majority spin 
band gap is not detached from the host valence band (cf. 
Refs. [13 ( 13 ) its pressure behavior is similar to that 
of the top of the valence band of the host material. 

Finally, in Fig. |2] we plot the pressure depen¬ 
dence of the majority spin band gap for GaAs and 



pressure (GPa) 

FIG. 2. (color online). Majority spin band gap of Gao .9375 
Mno,o 625 As and the band gap of GaAs plotted as a function 
of pressure. In both cases the direct-to-indirect band gap 
transition occurs at GPa. 


Gao.garsMng 0625 AS. GaAs has a direct band gap at am¬ 
bient conditions. Under hydrostatic pressure, the up¬ 
ward shift of the conduction band minimum at the F 
point eventually intersects the X minimum, and the ma¬ 
terial undergoes a direct-to-indirect band gap transition 
observed experimentally at ~4 GPa [13 • In Fig. we 
can see that the pressure behavior of the band gap of 
Gao, 9375 Mno.o 625 As is very similar to that of pure GaAs 
with a direct-to-indirect band gap transition at ~6 GPa. 


B. Cluster formation 

We consider the case of the formation of clusters com¬ 
posed of two Mn ions in GaAs. It is well established 
[SSHSg that the magnetic ions may occupy both substi¬ 
tutional cation sites, Mns, and As- and Ga-coordiiiated 
tetrahedral interstitial sites, Mnj(As) and Mnj(Qa), re¬ 
spectively. 

We start first with the configuration where the in¬ 
terstitial and substitutional ions are as far away from 
each other as possible within the supercell. For this 
case, we obtain that (Mns, Mnijca)) is more stable than 
(Mns,Mni(As)) by 0.11 eV. This result is congruent with 
previous computations 

Let us now focus on the formation of Mns“Mns and 
Mns-Mni complexes in GaAs. In both types of con¬ 
sidered clusters, the Mn ions are the nearest neigh¬ 
bors (NNs). The dissociation energies for Mns-Mns, 
Mns-Mni(As)i and Mns^Mnqca) pairs are defined as 

^L'ss = ii/(Mns,Mns) — ii/(Mns~Mns) 

ALIsi(As) = ■E(Mns,Mni(Ga)) - S(Mns-Mni(As)) (1) 

^^SI(Ga) = ^'(Mns, Mni(Ga)) — ^(Mas~Mni(Ga)) 

respectively, where E is the total energy of the 32- 
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TABLE I. Dissociation energies, interionic distances, and to¬ 
tal energy differences per ion between FM and AFM states 
for various configurations of Mn complexes in GaAs at zero 
pressure. Results for an isolated Mn dimer (Muj) are also 
included. 



dissociation energy 
(eV) 

d (A) 

Efm — Eafm 
(meV/Mn) 

Mns-Mni(Ga) 

0.42 

2.699 

69 

Mns-Mni(Aa) 

0.18 

2.879 

33 

Mn2 

0.21 

3.243 

30 

Mns-Mns 

0.11 

3.994 

-39 


or 33-atom supercell. To calculate iil(Mng, Mns) and 
iil(Mns, Mni(Ga))) we place the magnetic ions as far 
away from each other as possible within the super¬ 
cell. The largest possible distance between the mag¬ 
netic ions is 5.701 and 6.179 A for (Mng, Mns) and 
(Mns, Mni(Ga))j respectively. The calculated dissociation 
energies are 0.11, 0.18, and 0.42 eV for AEss, Ai?si(As)j 
and Aiilsi(Ga)) respectively. Our results suggest that 
pairs involving Mni are more likely to be present than 
purely substitutional clusters. It should be also noted 
that clusters involving Ga-coordinated interstitials are 
more stable than those involving As-coordinated inter¬ 
stitials. All mentioned values and some additional de¬ 
tails for the Mn complexes in GaAs are summarized in 
Table U where for comparison, we also include the iso¬ 
lated Mn 2 dimer for which we predict an antiferromag¬ 
netic (AFM) coupling with a binding energy of 0.21 eV in 
agreement with more sophisticated calculations involving 
quantum-chemical methods [3^ . The activation energy 
for the out diffusion of charged Mni ions that are bound 
in interstitial-substitutional complexes is a sum of the 
dissociation energy plus the energy barrier for diffusion 
of charged Mnj ions between interstitial sites [s^ . It 

should be, however, noted that the computed by us dis¬ 
sociation energies for clusters involving neutral Mni ions 
are equal to activation energies because of the lack of 
the migration barrier for swapping of neutral Mni ions 
between interstitial sites (3^ . 

Finally, in Fig. [3l we show the Mn-Mn NN distance d 
as a function of pressure. It is clear from the figure that 
the variation with pressure of d(Mns-Mni(Ga)) is much 
smaller than that of (i(Mns-Mns). For comparison, at 
zero pressure the unrelaxed values of (i(Mns-“Mns) and 
(i(Mns“Mni(Ga)) -where atoms are placed into ideal lat¬ 
tice positions- are 4.032 and 2.469 A, respectively (the 
relaxed values are listed in Table |I|. 


C. Curie temperature 

It is well established [B| that one of the factors limit¬ 
ing high-temperature ferromagnetism is the presence of 
Mns-Mni complexes in GaAs in which the Mn ions ex¬ 
hibit an AFM coupling. This coupling is not only pre- 



FIG. 3. (color online). Nearest neighbor distances d plot¬ 
ted as a function of pressure for Mns-Mni(Ga) (circles) and 
Mns-Mns (squares) pairs in GaAs. 


served but tends to increase in strength with pressure (for 
pressures up to ^7 GPa). This is shown in Fig. 01 where 
we plotted the total energy difference, E-pu — Eafm, be¬ 
tween the FM and AFM configurations for Mns-Mni(Ga) 
and Mng-Mng pairs as a function of pressure. For the 
substitutional Mn ions, we have considered NN and next 
nearest neighbor (NNN) positions in the Ga sublattice. 

The energy difference Sfm — Eafm for the substitu¬ 
tional ions can be used to evaluate the exchange inter¬ 
action Jij between i and j-site local spins. The ex¬ 
change energy for a system of interacting atomic mo¬ 
ments is given by the effective classical Heisenberg Hamil- 



pressure (GPa) 

FIG. 4. (color online). Total energy difference per ion be¬ 
tween the FM and AFM configurations for Mns-Mni(Ga) (cir¬ 
cles) and Mns-Mns (squares) pairs in GaAs as a function of 
pressure. For Mns-Mns, we have considered NN and NNN 
positions in the Ga sublattice. 
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pressure (GPa) 

FIG. 5. (color online). Curie temperature of (Ga,Mn)As cal¬ 
culated within MFA (squares) and MG (circles) methods as 
a function of pressure. The given values of the pressure co¬ 
efficients correspond to the region below 6 GPa marked by a 
dashed line. 

tonian H = where 

in our case Si = Sj = S = 5/2. For ferromag- 
netically coupled Mn ions at NN substitutional posi¬ 
tions —2 Jnn 5'^ = — ^'AFM- For zero pressure 

Efm — Eafm = —39 meV/Mn, therefore we obtain 
that Jnn'S'^ = 20 meV; the calculated in a similar way 
•/nnN'S'^ for NNNs gives a modest value of 5 meV [40| . 
The Curie temperature in the MFA is given by 

Tc = (2/3fcB)-(xE,z,Jo,5") = ^2) 

= (2S'^/3fcB) • (12 Jnn + SJnnn) /8 

where x is the concentration of Mn and zj denotes the 
number of sites on the j-th shell. The Curie temperature 
versus pressure is plotted in Fig. [5l From the slope of 
the MFA curve, we can estimate the derivative of tem¬ 
perature with respect to pressure dTc/dp = 7 K/GPa 
for pressures close to zero. A linear dependence of the 
MFA curve extends, however, up to ~6 GPa. Assum¬ 
ing a linear dependence of the Curie temperature ver¬ 
sus manganese concentration, we can estimate that for a 
6.25% concentration of Mn dTc/dp = 3.5 K/GPa. Be¬ 
yond 6 GPa the Curie temperature increases even faster, 
reaching the room temperature for ^7 GPa. 

We have also carried out MG simulations using the 
same exchange parameters Jnn and Jnnn as deter¬ 
mined above. From Fig. we can see that the pres¬ 
sure behavior of Tc calculated by means of MG is simi¬ 


lar to that obtained from MFA calculations. However, 
the critical temperature is shifted down, for all pres¬ 
sures, by an average value of 26%. For pressures close 
to 0 GPa, dTc/dp = 4.3 K/GPa , and an estimated 
value for x = 6.25% would be dTc/dp = 2.2 K/GPa, 
which is in the range of values obtained experimentally 
dTc/dp = (2-3) K/GPa at p < 1.2 GPa and estimated 
within the p-d Zener model for Gao.gsMno.orAs (l^ . 


IV. CONCLUSIONS 

In this paper, we have presented the pressure depen¬ 
dence of the Curie temperature for Gai_a;Mna;As with an 
experimentally realistic Mn contents of x = 6.25% and 
X = 12.5%, using the hybrid HSE06 functional for the 
description of exchange and correlation effects. In agree¬ 
ment with photoemission experiments, we have found for 
X = 6.25% and ambient pressure that the non-hybridized 
Mn-3d levels and Mn-induced states reside about 5 and 
0.4 eV below the Fermi energy, respectively. For ele¬ 
vated pressures, the Mn-3d levels, Mn-induced states, 
and the Fermi level shift towards higher energies, how¬ 
ever, the position of the Mn-induced states relative to 
the Fermi energy remains constant due to hybridization 
of the Mn-3c! levels with the valence As-4p orbitals. It 
follows from our computations that the Mn-derived spin- 
polarized feature in the majority spin band gap is not 
detached from the host valence band and its pressure be¬ 
havior is similar to that of the top of the valence band 
of the host material. We have also found that Tc at 
zero pressure is 181 and 258 K for MG and MFA cal¬ 
culations for X = 12.5%, respectively, while increases 
linearly under pressures up to 6 GPa. The estimated 
for X = 6.25% pressure-induced changes in Tc are +2.2 
and +3.5 K/GPa for MG and MFA calculations, respec¬ 
tively. The determined values of dTc/dp compare well 
with those found experimentally and estimated within 
the p-d Zener model for Gao.gaMno.orAs in the regime 
where hole localization effects are of smaller importance 
0 - 
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